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We study the implication of the recent measurement by the BELLE collaboration of 
the averaged branching fraction B exp [B -> (p,u) 7] = (1.8±g;| ± 0.1) x 10~ 6 for the CKM 
phenomenology. Combined with the averaged branching fraction B exp (B — > K* ^) = 
(4.06 ±0.26) x 10" 5 measured earlier, this yields R cxp [(p, u) 7/^*7] = (4.2± 1.3)% for the 
ratio of the two branching fractions. Updating earlier theoretical analysis of these decays 
based on the QCD factorization framework, and constraining the CKM-Wolfenstein pa- 
rameters from the unitarity fits, our results yield B t h[B — > (p, u) 7] = (1.38 ± 0.42) x 10~ 6 
and Rth[(p,uj) -y /K*j] = (3.3 ± 1.0)%, in agreement with the BELLE data. Leav- 
ing instead the CKM-Wolfenstein parameters free, our analysis gives (at 68% C.L.) 
0.16 < |V^/V^ S | < 0.29, which is in agreement with but less precise than the indirect 
CKM-unitarity fit of the same, 0.18 < | V^^/ V* s | < 0.22. The isospin-violating ratio in 
the B — > fry decays and the £?7(3)-violating ratio in the B® — > (p°,c<j)7 decays are pre- 
sented together with estimates of the direct and mixing-induced CP-asymmetries in the 
B — > (p, u) 7 decays within the SM. Their measurements will overconstrain the angle a 
of the CKM-unitarity triangle. 
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1. Introduction. Recently, the BELLE collaboration have presented evidence for the 
observation of the decays B + — > p + 7, B d — > p°7 and B d — > wy (and their charged 
conjugates) [1]. Their observation based on an integrated luminosity of 140 fb _1 lacks 
the statistical significance in the individual channels, but combining the data in the three 
decay modes and their charged conjugates yields a signal at 3.5a C.L. [1]: 

B exp [B - (p,u)i\ = (1.8±g;l ± 0.1) x 10~ 6 . (1) 

This result updates the previous upper bounds [2] by the BELLE collaboration, while the 
upper bound from the BABAR collaboration (at 90% C.L.) [3]: 

4xp[5^(p,w)7]<1.9xHr 6 , (2) 
remains to be updated. The experimental averages given above are defined as: 



B[B - (p, u>) j] = \ \B(B + -> p+ 7 ) + T -f- [B(B° d -> A) + B(B° d -> ury)] } 



(3) 



and the world average [4] for the 5-meson lifetime ratio: 

t b+ /t b o = 1.086 ±0.017, (4) 

has been used in arriving at the BELLE result This is the first observation of 

the CKM-suppressed electromagnetic penguin b —>■ d'j transition. The CKM-allowed 
b — > S7 transition in the exclusive decays B — > K*j was observed more than a decade 
ago by the CLEO collaboration [5], followed by the observation of the inclusive decay 
B — > A s 7 in 1994 [6]. Since then, data on these decay modes have been provided by 
a number of experimental collaborations, and the current situation is summarized in 
Tabled In getting the isospin-averaged branching ratio B exp (B — ► K*^), we used the 
following definition: 



B cxp (B - iT 7 ) = ~ 



B cxp (B + - JC+ 7 ) + ^ B cxp (B° d - K*°y) 



(5) 



and the world average (j3J) for the 5-meson lifetime ratio. Table ^ also contains the 
measurements of the inclusive decay B — > X s 7 branching fraction, the resulting ratio 
of the exclusive-to-inclusive decay rates R exp (K*^ / X s j), for each experiment separately, 
and their world averages, with the errors added in quadrature. 

The measurements from BELLE and the upper limit from BABAR on the B — > 
(p, uj)^ decays given in (JTJ) and (J2J), respectively, can be combined with their respective 
measurements of the B — > K*^ decay rates to yield the following ratios: 

R cxp [(p, uh/K^} < 0.047, (BABAR) (6) 

R«xp[(p,w)'y/K*i\ = 0.042 ±0.013, (BELLE) (7) 

where R exp [(p, u)^/K*^} = B cxp [B — ► (p, cj) 7]/S oxp (_B — > K*j). In this paper, we do an 
analysis of the two quantities in Eqs. (0) and (J7|) in the context of the SM. 
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Table 1: The branching ratios averaged over the charge-conjugated modes (in units 
of 10~ 5 ) of the exclusive decays B + — > K* +, ~f and B d —>■ K°* / y and the inclusive decay 
B — > X s 7 taken from Refs. [4,7-10]. The averaged branching ratios defined in © and 
the ratio of the exclusive-to-inclusive branching ratios R exp {K* r y / 'X s 7) are also tabulated. 



Fraction 


BABAR 


BELLE 


CLEO 


Average 


B cxp (B+ - K*+ 7 ) 
B exp (B° d - *T*° 7 ) 
B cxp (B K* 7 ) 


3.87 ±0.28 ±0.26 
3.92 ±0.20 ±0.24 
4.06 ± 0.26 


4.25 ±0.31 ±0.24 
4.01 ±0.21 ±0.17 
4.30 ±0.25 


3.76j:gJ| ± 0.28 
4.55i^± .34 
4.35 ± 0.62 


4.03 i 0.26 
4.01 i 0.20 
4.20 i 0.17 


£exp(# - X s7 ) 
R exp (K*j/X sl ) 


38.8 ± 3.6±| | 
105+ ' 021 


35.5 ± 3.2±|;S±5;^ 

121+ 0019 
u.j.zi_ 015 


32.li4.3jJJ 

n 1 o fi +0.033 


35.1 i 3.0 
0.117 i 0.012 



2. Effective Hamiltonian. The starting point for the theoretical discussion of the 
radiative b — > decays (equivalently B ^ fry and B — > u^y decays) is an effective 
Hamiltonian obtained from the Standard Model (SM) by integrating out the heavy degrees 
of freedom (the top quark and W /± -bosons). The resulting expression at the scale /i = 
0(m&), where mi, is the 6-quark mass, is given by 



H 



eff 



G f 

71 



V ubV* d 



+ v cb v: d 



ct\n)of\n) 



- v tb v; d [cf(/i) om + cf(n) om] +...}, 

where Gp is the Fermi coupling constant, and only the dominant terms are shown. The 
operators and , (q — u,c), are the standard four-fermion operators: 



0[ q) = (d altl (l- l5 )qp) (q^(l-J 5 )b c 



Of 



(d a ^(l-J 5 )q a ) (^(1-75)^), (9) 



and Oj and 0$ are the electromagnetic and chromomagnetic penguin operators, respec- 
tively: 



o 7 



errib 



(d a a^(l + l5 )b a )F f 



Os = |^ (4^(1 + 75) W G%. (10) 



Here, e and g s are the electric and colour charges, F^ u and G a ^ v are the electromagnetic and 
gluonic field strength tensors, respectively, T?g are the colour SU (N c ) group generators, 
and the quark colour indices a and f3 and gluonic colour index a are written explicitly. 
Note that in the operators O-? and Og the <i-quark mass contributions are negligible and 



therefore omitted. The coefficients C[ q \fi) and C^\n) in Eq. (jHJ) are the usual Wilson 



coefficients corresponding to the operators O^ and Of , while the coefficients Cj S ([i) 



>(</) 
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and C| ff (p) include also the effects of the QCD penguin four-fermion operators which 
are assumed to be present in the effective Hamiltonian (JSJ) and denoted by ellipses there. 
For details and numerical values of these coefficients, see [11] and reference therein. We 
use the standard Bjorken-Drell convention [12] for the metric and the Dirac matrices; 
in particular 75 = ij ^ 1 ^ 2 ^ 3 , and the totally antisymmetric Levi-Civita tensor e^ upa is 
defined as e i23 = +1- 

For the b — > 37 decay (equivalently the B — > K*^ decays), the effective Hamilto- 
nian 'H b e ^ s describing the b — > s transition can be obtained by the replacement of the 
quark field d a by s a in all the operators in Eqs. (JHJ) and (fTU|) and by replacing the CKM 
factors V qb V* d — > V qb V* s (q = u,c,t) in H^ d ©• Noting that among the three fac- 
tors V q bV* , the combination V ub V* s is CKM suppressed, the corresponding contributions 
to the decay amplitude can be safely neglected. Thus, within this approximation, unitar- 
ity of the CKM matrix yields V c bV* s = —V t bV t * s , the dependence on the CKM factors in the 
effective Hamiltonian Ti^ s factorizes, and the CKM factor is taken as VtbV t *. Note also 
that the three CKM factors shown in l~L h e ^ d are of the same order of magnitude and, hence, 
the matrix elements in the decays b — > d'j and B — > (p, 07)7 have non-trivial dependence 
on the CKM parameters. 

3. Theoretical framework for the B — ► decays. To get the matrix elements for 
the B — > V'-y (V = K*,p,u) decays, we need to calculate the matrix elements (Vj\Oi\B), 
where Oi are the operators appearing in 'H b ^ s and 7~i b e ^ d . At the leading order in a s , 
this involves only the operators O7, and where the latter two are important 

only for the B — > (p, 07)7 decays. One also uses the terminology of the short- distance 
and long-distance contributions, where the former characterizes the top-quark induced 
penguin-amplitude and the latter includes the penguin amplitude from the u- and c-quark 
intermediate states and also the so-called weak annihilation and VT-exchange contribu- 
tions. There are also other topologies, such as the annihilation penguin diagrams, which, 
however, are small. For a recent discussion of the long-distance effects in B — > V'-f decays 
and references to earlier papers, see Ref. [13]. 

Including the 0(a s ) corrections, all the operators listed in (JSJ) and (fTU|) have to be 
included. A convenient framework to carry out these calculations is the QCD factorization 
framework [14] which allows to express the hadronic matrix elements in the schematic 
form: 

1 

<yi\Oi\B) =F B - v T i I + l^±J du^k^ik+^Ku), (11) 



where F B ^ V are the transition form factors defined through the matrix elements of the 
operator C 7 , B + (/c + ) is the leading-twist 5-meson wave-function with k + being a light- 
cone component of the spectator quark momentum, <fi\{u) is the leading-twist light-cone 
distribution amplitude (LCDA) of the transversely-polarized vector meson V, and u is 
the fractional momentum of the vector meson carried by one of the two partons. The 
quantities % ! and T* 1 are the hard-perturbative kernels calculated to order a s , with 
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the latter containing the so-called hard-spectator contributions. The factorization for- 
mula (|11|) holds in the heavy quark limit, i.e., to order Aqcd/^b- This factorization 
framework has been used to calculate the branching fractions and related quantities for 
the decays B — > K*^ [15-17] and B —>■ p7 [15,17]. The isospin violation in the B — > K*^ 
decays in this framework have also been studied [18]. (For applications to B — > K*'-/*, see 
Refs. [16,19,20]). Very recently, the hard-spectator contribution arising from the chromo- 
magnetic operator 8 have also been calculated in next-to-next-to-leading order (NNLO) 
in a s showing that the spectator interactions factorize in the heavy quark limit [21]. How- 
ever, the numerical effect of the resummed NNLO contributions is marginal and we shall 
not include this in our update. 

In what follows we shall use the notations and results from Ref. [15], to which we 
refer for detailed derivations, and point out the changes (and corrections) that we have 
incorporated in this analysis. The branching ratio of the B — > K*^ decay corrected to 
0(a s ) can be written as follows [15]: 

(12) 

where a is the fine-structure constant, m 6 pole is the 6-quark pole mass, and M B and m K * 
are the B- and i^*-meson masses, respectively. The quantity ^ ^ is the soft part of 
the QCD form factor T-^* (q 2 ) in the B — > K* transition, which is evaluated at q 2 = 
in the HQET limit. For this study, we consider ^ ' as a free parameter; its value will 
be extracted from the current experimental data on B — > K*^ decays. Note that the 
quantity ^ ^ used here is normalized at the scale \x = m^ po i e of the pole 6-quark mass. 
The corresponding quantity in Ref. [19] is defined at the scale \i = m^pg involving the 
potential-subtracted (PS) 6-quark mass [22,23], which is numerically very close to the 
pole mass used here. 

The function C^ efi (fi) in Eq. ()12j) is the Wilson coefficient of the electromagnetic 
operator Oj in the leading order and the function A^> (/x) includes all the NLO corrections: 

= <V) + A£J(A0 + AgT^p), (13) 

where Aq . A^ T and A^ K denote the 0(a s ) corrections in the Wilson coefficient Cf 5 , 
the b — > sj vertex, and the hard-spectator contributions, respectively. Their explicit ex- 
pressions are given in Eqs. (5.9), (5.10) and (5.11) of Ref. [15]. The values used in the 
numerical analysis are collected in Table El Some comments on the input values are in 
order. The top-quark mass (interpreted here as the pole mass) has been recently updated 
and revised upwards by the Tevatron electroweak group [24], and the new world aver- 
age m t ,po\e = (178 ± 4.3) GeV is being used in our analysis. The product IV^V^I of the 
CKM matrix elements can be obtained from the estimate |T^&| = 0.0412 ± 0.0021 [25] 
using the relation \V tb V*\ ~ (1 - X 2 /2)\V cb \, which yields \V tb V*\ = 0.0402 ± 0.0020 for 
A = 0.2224. The S"?7(3)-breaking effects in the K- and X*-meson LCDAs have been 
recently re-estimated by Ball and Boglione [26]. In this update, the transverse decay 



B th {B - iT 7 ) = r B 



v th v t 



ts I 



32 ii 4 



;{«*) 



m K* 

Ml 



5 



Table 2: Input quantities and their values used in the theoretical analysis. The values of 
the masses, coupling constants and given in the first four rows are fixed, and those of 
the others are varied in their indicated ranges to estimate theoretical uncertainties on the 
various observables discussed in the text. 



Parameter Value 


Parameter Value 


M w 80.423 GeV 

M B 5.279 GeV 

G F 1.16639 x 10" 5 GeV~ 2 

a s {M z ) 0.1172 


M z 91.1876 GeV 
m K * 894 MeV 
a 1/137.036 
A h 0.5 GeV 


m t , po ie (178.0 ± 4.3) GeV 
\V tb V*\ (40.2 ± 2.0) x 10- 3 
f B (200 ± 20) MeV 
afP (1 GeV) -0.34 ±0.18 
\ B ] + (1 GeV) (2.15 ± 0.50) GeV' 1 


m b , po ie (4.65 ± 0.10) GeV 
y fi = m c /m b 0.27 ±0.06 

/[^ (1 GeV) (182 ± 10) MeV 
af* ] (1 GeV) 0.13 ±0.08 
a B ,+ (1 GeV) 1.4 ±0.4 



constant of the i^*-meson, /j^ \ has remained practically unchanged, but the Gegen- 
bauer coefficients in the i^*-meson leading-twist LCDA are effected significantly. The 
two Gegenbauer moments ■* and ^ used in the calculation of the hard-spectator 
contributions are now larger in magnitude, have larger errors and, moreover, the first 
Gegenbauer moment changes its sign. For comparison, previously, these coefficients were 
estimated as af*\l GeV) = 0.20 ± 0.05 and af*\l GeV) = 0.04 ± 0.04. The effect of 
these modifications on the QCD form factor T^*(0), as well as of some other technical 
improvements [26], has not yet been worked out. Lastly, the first inverse moment \ : B 1 + (fi) 
of the B- meson LCDA has also changed. In our previous analysis [15], we used the value 
A5+(/i sp ) = (3.0 ± 1.0) GeV- 1 where the error effectively includes the scale dependence 
of the leading-twist light-cone 5-meson wave-function <pB,+{k, ft). In a recent paper by 
Braun et al. [27], the scale dependence of this moment is worked out in the NLO with the 
result: 

A B » = A^> ) (l - L, + (AK>) - l\ I*-) , (14) 

I n L 1 J ) 

where (a s Cp/T!') ln^/zio) < 1 and the quantities A^ 1 + ( / u) and a B + {n) are defined as 
follows: 



oo 



{^) = J ^T<pB,+( k ^)' = A B)+ (//) J y ln^0 Bi+ (fc,/x). (15) 







At the initial scale fio = 1 GeV of the evolution, the above quantities were estimated by 
using the method of the Light-Cone-Sum-Rules (LCSR) and their values are presented 
in Table El At the typical scale fi sp = A/Ah^poic — 1-52 GeV (here, A^ = 0.5 GeV is 
a typical hadronic scale) of the hard-spectator corrections, the first inverse moment is 
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now estimated as: A^ 1 + (/i sp ) = (2.04 ±0.48) GeV -1 . Note that, while overlapping within 
errors with the previously used value, the updated estimate is substantially smaller as 
well as the current error on this quantity is now reduced by a factor of two. 

Updating the analysis presented in Ref. [15], and using the experimental results on 
the branching ratios for the B — * K*y and B — > X s y decays given in Table ^ the 
phenomenological values of the soft part of the QCD form factor are: ^ \o) = 0.28 ± 

0.02, ^* ±} (0) = 0.27 ±0.02 and £ { ** /Xs \o) = 0.25 ±0.02 resulting from the B° d -> K*°y 
and B ± — > K ^ branching ratios and from the ratio R exp (K*'y / Xgj), respectively. The 
QCD form factor Tf*(0) differs from its soft part ^ \o) by 0(a s ) terms worked out in 
Ref. [19], which in our notation is given in Eq. (5.13) of Ref. [15]. However, the updated 
input parameters reduce this correction, yielding typically a correction of 2 — 4% only, 
in contrast to about 8% previously. Thus, the QCD transition form factor (0) now 
differs only marginally from its soft part, and is estimated as follows: 

ff(0) = 0.27 ±0.02. (16) 

The central value of the QCD form factor (|16|) extracted from the current data has re- 
mained unchanged compared to the previous estimate (0) = 0.27±0.04 (see Eq. (5.25) 
of Ref. [15]), but the error is now reduced by a factor 2, mostly due to the reduction of 
the uncertainty on the input parameters. It remains an interesting and open theoretical 
question if improved theoretical techniques for the calculation of the transition form factor 
T/ £: *(0) could accommodate this phenomenological result. 



4. Results for B — > (p, uj) 7 decays and comparison with the BELLE data. This 
part is devoted to an update of the theoretical predictions for the B — > /ry and B® w-f 
branching ratios, and their comparison with the BELLE data. Results for the direct and 
mixing-induced CP-violating asymmetries in these decays, the isospin-violating ratio in 
the B — > p7 decays, and the S'[/(3)-violating ratio in the neutral B® — > p°7 and B% — > 
decays are also presented. 



4.1. Branching ratios. We now proceed to calculate numerically the branching 
ratios for the — > p ± 7, B® — > p°7 and B d decays. The theoretical ratios 

involving the decay widths on the r.h.s. of these equations can be written in the form: 



R th { P l/K^) 



B th (B -> pry) 



B th (B 



V* 



td 



(M 



m 



2N3 



V ts 

v td 



Is 



(M 



m 



2 \3 



(M- 



m 



K 



: 



e [1+AR( P /K*)\, (17) 
C 2 [l + AR(u/K*)}, (18) 



where m p and m w are the masses of the p- and cj-mesons, ( is the ratio of the transition 
form factors, £ = Tf(0)/T 1 K *(0), which we have assumed to be the same for the p°- and u- 
mesons, and S p — 1 and 1/2 for the p ± - and p°- meson, respectively. To get the theoretical 
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Table 3: Input parameters and their values used to calculate the branching fractions in 
the B — > p7 and B® — ► decays. The parameters entering in the B —>■ K*^ part in 
Eqs. (JT7J) and (JTSJl are given in Tabled 



Parameter Value 


Parameter Value 


m p 771.1 MeV 


m w 782.57 MeV 


GeV) (160 ±10) MeV 
C 0.85 ±0.10 
4 ±} ±0.30 ± 0.07 
p 0.17 ±0.07 


a ( H(l GeV) 0.20 ±0.10 
|V^* d | (8.1 ±0.8) x 10- 3 
ej) = - e ^ } +0.03 ± 0.01 
fj 0.36 ±0.04 



branching ratios for the decays B — > pj and B® — > uj, the ratios ()17|> and (j!8)l should be 
multiplied with the corresponding experimental branching ratio of the B —>■ K*^ decay. 

The theoretical uncertainty in the evaluation of the -Rth(p7 / 'K*j) and R t h(uJj/K*j) 
ratios is dominated by the imprecise knowledge of ( = T 1 p (0)/T/ s '*(0) characterizing the 
SU(3) breaking effects in the QCD transition form factors. In the S"t/(3)-symmetry limit, 
^f(0) = ^i *(0)> yielding £ = 1. The SU (3)-bieakmg effects in these form factors have 
been evaluated within several approaches, including the LCSR and Lattice QCD. In the 
earlier calculations of the ratios [15,28], the following ranges were used: ( = 0.76±0.06 [15] 
and ( = 0.76 ±0.10 [28], based on the LCSR approach [29-33] which indicate substantial 
SU(3) breaking in the B — > K* form factors. There also exists an improved Lattice 
estimate of this quantity, ( = 0.9±0.1 [34]. In the present analysis, we use ( = 0.85±0.10, 
given in Table El together with the values of the other input parameters entering in the 
calculation of the B — > (p, to) 7 decay amplitudes. 

We now discuss the difference in the hadronic parameters involving the p°- and 00- 
mesons. It is known that both mesons are the maximally mixed superpositions of the uu 
and dd quark states: |p°) = (\dd) — \uu))/y/2 and \uo) = (\dd) + \uu))/y/2. Neglecting the 
^-exchange contributions in the decays, the radiative decay widths are determined by 
the penguin amplitudes which involve only the \dd) components of these mesons, leading 
to identical branching ratios (modulo a tiny phase space difference). The M^-exchange 
diagrams from the 0^ and operators (in our approach, we are systematically 

neglecting the contributions from the penguin operators O3, yield contributions 

equal in magnitude but opposite in signs. In the numerical analysis, the LCSR results: 

= ±0.03 ± 0.01 and = —0.03 ± 0.01 [30], are used, where the smallness of these 
numbers reflects both the colour-suppressed nature of the W-exchange amplitudes in 
B% — > (p°,u;)7 decays, and the observation that the leading contributions in the weak 
annihilation and W -exchange amplitudes arise from the radiation off the <i-quark in the 
-B°-meson, which is suppressed due to the electric charge. The parameter entering in 
B ± — > p ± 7 and ef in the B° d — > p°7 decay have been estimated in the factorization ap- 
proximation for the weak annihilation (and W-excliange) contribution, but this is expected 
to be a good approximation in the heavy quark limit, where the 0(a s ) non-factorizable 
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Table 4: Updated theoretical estimates of the functions AR(p/K*) and AR(u/K*), and 
the ratios of the branching ratios R^p^y/K*^) and Rt^wy / K*j) defined in Eqs. ()17|) 
and (JUJ), respectively. The third and fourth rows give the branching ratios B t h(B — > /ry) 
and £> t h(-B° — > ory) (in units of 10 -6 ) and direct CP asymmetries in the B — > p7 and 
5° — >■ CU7 decays, respectively. 





fi ± -> p± 7 


B° d - A 


B° d - u; 7 


Ai? 


0.116 ±0.099 


0.093 ±0.073 


0.092 ±0.073 




0.0334 ±0.0103 


0.0164 ±0.0049 


0.0163 ±0.0049 


Bo, 


1.35 ±0.42 


0.66 ±0.20 


0.65 ±0.20 


/Idir 


(-11.6 ±3.3)% 


(-9.4+^)% 


(-8.811*)% 



corrections are found to be suppressed in the chiral limit [13]. Moreover, their magnitudes 
can be checked experimentally through the radiative decays B^ — > £1^7, as emphasized 
in Ref. [13]. These and the other parameters needed for calculating the branching ratios 
in the B — > (p, u) 7 decays are given in Table 01 where we have also given the default 
ranges for |VtbV^| and the CKM-Wolfenstein parameters p and fj obtained from the recent 
fit of the CKM unitarity triangle [25]. 

The individual branching ratios B t h(B — > pry) and B t h(B^ — > ujy) and their ratios 
Rth{pj/ K*j) and Rt^usy / K*y) with respect to the corresponding B —> K*y branching 
ratios are presented in Table El Note that in our estimates there is practically no difference 
between the B® — > p°7 and B® — > 077 branching fractions, as the two differ only in the signs 
of the weak-annihilation contributions in the decay amplitudes, but these contributions 
given in terms of the parameters and e£ are small. Using the definition of the 
weighted average (jSJ), we get: 

B th [B -> (p,cj) 7 ] = (1.38 ± 0.42) x 10" 6 , (19) 
R\ h [(p,u)j/K*j) = 0.033 ±0.010, (20) 

where the current experimental values of the B — > K* r y branching ratios given in Table ^ 
have been used in arriving at the result (|19|). These theoretical estimates, carried out 
in the context of the SM, are in the comfortable agreement with the current BELLE 
measurements (JTJ) and (J7J). 



4.2. CP-violating asymmetries. The direct CP-violating asymmetries in the 
decay rates for B + — > p +r y and B® — > (p°,u)'y decays and their charged conjugates are 
defined as follows: 

.dw + x _ B{B~ - p- 7 ) - B{B+ - p+ 7 ) 



B[B- -> p-j) + B{B+ -> p+7) ' 



A^(pS) = ^ A) = B{B ° - A) (21) 
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The explicit expressions for the first two of these asymmetries in terms of the individual 
contributions in the decay amplitude can be found in Ref. [15] and the one for the last, 
Aqp(ujj) may be obtained from A^>{p j) by obvious replacements. Their updated values 
in the SM, taking into account the parametric uncertainties and adding the various errors 
in quadrature, are presented in Table EJ The main contribution to the errors is coming 
through the scale dependence and the uncertainty in the c- to 6-quark mass ratio, which 
is a NNLO effect. A complete NNLO calculation will certainly be required to reduce the 
theoretical errors. It should be noted that the predicted direct CP-asymmetries in all 
three cases are rather sizable (of order 10%) and negative. This differs from our earlier 
estimates [15,28], worked out for ^4cp(p ± 7) an d *4cp(p°7)> where the explicit expressions 
were erroneously typed and used in the numerical program with the incorrect overall sign. 

The dependence of the direct CP-asymmetry on the CKM unitarity-triangle angle a 
is presented in the left frame in Fig. ^ We note that the CP-asymmetries are calculated 
with the strong phases generated perturbatively in 0(a s ) in the QCD factorization ap- 
proach. In particular, they do not include any non-perturbative rescattering contribution. 
We recall that for the CP- asymmetries in non-leptonic decays, such as in B — > tttt, current 
data point to the inadequacy of the perturbatively generated strong phases [25]. In radia- 
tive decays B — > (p, u)^, such long-distance effects enter via the penguin amplitudes Pu\ 
which are the tm-loop contributions involving the operators 0\ (i = 1,2), and Pc \ the 
corresponding cc-loop contributions involving the operators [13]. They are included 
in the estimates of the complete matrix elements to a given order [here, up to 0(a s )]. 
In the hadronic language, they can be modelled via the hadronic intermediate states, 
such as 5 ± — > p^p — > p ± 7, B ± — ► £)*±i)* _ > p ± 7, etc. Their relative contribution 
at the amplitude level was estimated for the decay B~ — > p~7 as \P C /P t \ ~ 0.06 [13], 
with \P U \ <C \P C \. A recent model-dependent estimate [35] of the long-distance contri- 
bution in B° — > p°7 via the intermediate D + D~ state, B° — > D + D~ — > p°7, puts the 
relative contribution of the long-distance (LD) and short-distance (SD) contributions to 
the decay widths as Pld/Psd — 0.3, using the lowest order result for r SD . Taking into 
account that the next-to- leading order contributions in r SD , updated in this paper, result 
in an enhancement by a factor of about 1.7, and noting further that the perturbative 
charm-penguin contribution should be subtracted from Tld to avoid double counting, 
the remaining rescattering contributions are very likely below 10%. However, one can 
not exclude an enhanced charm-penguin contribution at this rate and the CP-asymmetry 
^4qp(p°t) cou ld beinfluenced from such long-distance contribution. Charm-penguin en- 
hanced effects can be also tested in the Dalitz pair reaction B° — > p 0/ j* — > p° e + e~ through 
measurements of the Stoke's vector components [35]. 

We now discuss the time-dependent (or mixing-induced) CP-asymmetry in the B®(t) — > 
(p°, uj) 7 and B®{t) — > (p°, uj) 7 decays. Below, the equations for the S^-meson decays into 
the final state with the p°-meson are presented. Similar quantities for the decays with the 
w-meson production can be obtained by the obvious replacement: — > ■ 
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Figure 1: Left frame: The direct CP- violating asymmetries for the B — > p7 and — > oo^y 
decays, defined in Eq. ()21j). Right frame: The mixing-induced CP-violating asymmetries 
for the B^ — > p°7 and B® —>■ decays, defined in Eq. (|22|) . 



The time-dependent CP-asymmetry in the decays of neutral 5°-mesons and its CP- 
conjugate involves the interference of the B^ — B® mixing and decay amplitudes and is 
given by [36]: 

G gp(t) = -C n cos(AM d t) + S n sin(AM d t), (22) 

where AM^ ~ 0.503 ps -1 is the mass difference between the two mass eigenstates in the 
B® — Bj system. For getting explicit formulae for C P1 and Spy, it is convenient to introduce 
the quantity: 

_ g A(B° d -> A) _ cf cff + - [cf cff ^ 0) + A»] Fe^ 

" P - P°7) cf eff + AW* - [cf + A«] F e- ' 1 } 

where p/q ~ exp(2i/3) is the 5° — B^ mixing parameter and F = Rb/Rt with R/, = 
\Jp 2 + ff and R t = a/(1 — p) 2 + ff- In terms of A P7 , the direct and mixing- induced 
CP-violating asymmetries can be written as follows: 

C n = -A&(p°y) = 1 —^ 2 , Spy = (24) 

1 -+- \Apy\ 1 T~ \Apy\ 

Thus, the direct CP-violating asymmetry C P1 is expressed by Eq. (6.6) in Ref. [15] while 
the mixing- induced CP-violating asymmetry Spy in NLO can be presented in the form: 

cNlo = oLO _ 2Fsina [1 - 2Fef cos a + (Fef) 2 cos(2a)] A% - e^A™ 

'7 



[l-2Fefcosa+(Fef) 2 } 2 ^ 0)eff 
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LO _ 2Fef sin a (1 - Fef cos a) 
1 — 2Fe\' cos a + (Fe A ') 2 

where A^ 1 and are the real parts of the NLO contributions to the decay amplitudes 
entering Eq. (f2*3*J) . It is easy to see that, neglecting the weak-annihilation contribution 
(e^ = 0), the mixing- induced CP-asymmetry vanishes in the leading order. However, 
including the 0(a s ) contribution, this CP-asymmetry is non-zero. 

The dependence on the CKM unitarity-triangle angle a of the mixing- induced CP- 
asymmetry for the B^-meson modes considered is presented in Fig. ^ (right frame). The 
dashed lines show the dependence in the LO while the solid lines correspond to the NLO 
result. Thus, fixing the parameters to their central values, one notices a marked effect from 
the NLO corrections on both and S^®. However, including the errors in the input 
parameters, the resulting allowed values for and are rather uncertain. This is 
worked out by taking into account the SM range a = (92 ± 11)° [25], and the numerical 
values for these asymmetries in the leading order and including the 0(a s ) corrections are 
as follows: 

5 p L 7 ° = (-2.7 ± 1.0)%, 5 p N 7 LO = (0.lil; 3 7 )%, (27) 
S™ = (+2.7 ± 1.0)%, = (4.011;°)%. (28) 

Thus, the ±lcr ranges for the mixing- induced asymmetries in the SM are: —0.04 < 
< 0.05 and —0.01 < < 0.09. They are too small to be measured in the near 

future. Hence, the observation of a significant (and hence measurable) mixing-induced 
CP-asymmetries S n and would signal the existence of CP-violating phases beyond 
the SM. 

4.3. Isospin-violating ratio. The charge-conjugated isospin-violating ratio is de- 
fined as follows: 

A = - [A +0 + A" ! , A ±0 = P±7) l. (29) 

The explicit NLO expression in terms of the vertex, hard-spectator and weak-annihilation 
contributions to the decay amplitude can be found in Ref. [15]. The dependence of this 
ratio on the angle a is shown in the left frame in Fig. El With the improved input, the 
updated result is: 

A = (1.1 ±3.9)%. (30) 

Thus, the isospin violation in B — > /ry decays is expected to be small in the SM. The 
reason for this lies in the dependence A oc cos a + 0[(e^) 2 , a s ], and we have used 
the current knowledge of the angle a from the CP-asymmetry in B — > titt decays and the 
indirect unitarity fits, yielding a = (92 ± 11)° [25]. 



12 




0.10 



0.08 
^ 0.06 

3 o.c 

0.02H 



0.00 



1 

- 


qa/t 

bM 


i 

' / ' 
' / ' 
' / ' 
' / ' 
' / ' 
if' 

i / ' 
' / ' 
' / ' 
/ / ' 
' / ' 


Rexp (BELLE) 


/ 

/ 

/ 


/ / 
/ / 
/ ' 




i / 


/ 


/ / 
/ / 
/ / ' 
' / ' 

J*' 

1 1 


J 





0.1 0.2 0.3 

\V U /Vt s \ 



0.4 



Figure 2: Left frame: The isospin-violating ratio A, defined in Eq. (|29j) for the de- 
cays B — > p7, and the S't/(3)-breaking ratio A^ p ^\ defined in Eq. (13 1|) and involving 
the B% — > (p°,c<j)7 decays, plotted as functions of the CKM unitarity-triangle angle a. 
Right frame: The ratio Rth[{p, uS)'y/K*'y\ plotted as a function of |Kd/Ks|- The current 
experimental measurement with its ±lcr range is shown as the horizontal band. The solid 
and dashed curves are the theoretical predictions in the SM and their ±lcr errors, respec- 
tively. The vertical dotted lines show the SM-based best-fit interval for |V^/Vi S | from the 
CKM unitarity fits. 



4.4. S77(3)-violating ratio. Another quantity of experimental interest is the ratio 
A&M^ involving the B® — ► (p, u)j decays. It can be defined as: 



with 



A 



A 



B "t" B 



(31) 



W») - 


(Ml 


-ml 


) 3 B(B° d 


- A) - (Ml 


-m 2 p 


) 3 B(B°- 


-> urf) 


B — 


(Ml 


-ml 


) 3 W 


- p° 7 ) + (Ml 


-m% 


3 B(B° d - 


-> uj) 


(pM - 


(Ml 


-ml 


fB(B° 


- A) - (Ml 


-m 2 p 




-> urf) 


B ~ 


(Ml 


-ml 


?B{B« 


- A) + (Ml 


-m 2 p ] 


3 B(B° d - 


-> uj) 



The weighted factors in A^^ are introduced to suppress the effect of the phase space 
due to the difference in the p- and cj-meson masses. The expression for A^^ , derived in 



13 



0(a s ) order, is rather lengthy: 



B,B 



cos a - FA U R =f Af ]t sin a 



Cf )cS [cosa-Fe A } 



(32) 



A { R ]t - FA U R cos a q= FAy sin a 



Cf )eff [1 - 2Fe A cos a + (Fe A ) 2 + (FAea) 2 ] 



where Sa = (^a + £ a')/^ an< ^ ^ £ a — ( £ a' ~ £ a')/^- our approximation, sa - 
neglecting tiny corrections ~ (FAsa) 2 , the final expression is greatly simplified: 



» 



-2FAe A [cos a - Fe A ] 



1 - 2Fe A cosa + (Fe^) 2 + (FAe A ) 



(o) 



A 



(pAO 

NLO 



2FAe, 



a 



(0)eff 



(Cf )eff - A^ 1 ) cosa ± FA^ cos(2a) 



i(i)t 



(33) 
and, 

(34) 



The dependence of this ratio on the angle a is shown in the left frame in Fig. |21 In the 
SM, with the input parameters specified above, this ratio can be estimated as: 



A 



(pM 

NLO 



(0.3 ±3.9) x 10" 



(35) 



This value is an order of magnitude smaller than the isospin- violating ratio ()30|) in B — > p7 
decays due to the suppression of the weak-annihilation contributions in the decays of the 
neutral F-meson. In this case, the neglected subdominant long-distance contributions 
may become important. They can be estimated in a mo del- dependent way. In any 
case, the result in (|33j) should be improved by including the contributions of the penguin 
operators and the NNLO corrections. The ratio A^ P//<JJ - ) in the SM is also too small to be 
measured. Both the ratios A and A^' w ' are sensitive tests of the SM, and as argued in 
Refs. [28,37] for the isospin-violating ratio A, their measurements significantly different 
from zero would reveal physics beyond the SM. 



5. Determination of |V^/V^| from R exp [(p, u) 7/^*7]. To extract the value of 
iVtd/VtsI from the B — > (K* , p,u) 7 decays, we use the ratio R^Kp, u) 7/^*7], which 
can be rewritten within the SM as follows: 



R th [(p,u) 1 /K* 1 ]=r 



(pM 

th 



c 2 



>/") 

th 



1.18 ±0.10, 



(36) 



where the error in r^/^ takes into account all the parametric uncertainties except in ( 
and |Vtd/Vt s | which are treated as free variables. Applying this equation to the BABAR 
upper limit (0) and the BELLE experimental range (j7J), the product C|Kd/Ks| can be 
restricted as follows: 



C\v td /v ts \ 

(\Vtd/V ts \ 



> 0.19, 

= 0.19 ±0.03, 



(BABAR) 
(BELLE) 



(37) 
(38) 
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where the bound is at 90% C.L., following from the BABAR data. At present, the error 
in (|38|) is dominated by the experimental uncertainty. Using the range ( = 0.85 ±0.10 for 
the ratio of the transition form factors, one gets the following constraints on the CKM 
matrix element ratio |Vjd/\4 s |: 

\Vtd/V u \ > 0.19, (BABAR) (39) 

\V td /V ts \ = 0.22 ± 0.05, (BELLE) (40) 

where the lower limit from the BABAR data © corresponds to 90% CL. In arriving at 
these numbers, the theoretical and experimental errors were considered as uncorrelated. 
Taking this correlation into account, the BELLE data yields the range 0.16 < jV^/V^I < 
0.29, which is much larger than but in agreement with the SM range |Vtd/^ts| = 0.20±0.02. 

The dependence of the ratio R^Hp, u)j / K*j] on |Vfd/Vt s | is shown in the right frame 
in Fig. El The solid curve corresponds to the central values of the input parameters, and 
the dashed curves are obtained by taking into account the ±lcr errors on the individual 
input parameters in R t h[{p, w)'j/K*'-f] and adding the errors in quadrature. The current 
measurement for this quantity is also shown in this figure. Experimental error is currently 
large which renders the determination of |V^/V^| uncertain. However, in the long run, 
with greatly increased statistics, the impact of the measurement of R exv [(p, lj)j/K*j] on 
the CKM phenomenology, in particular the profile of the unitarity triangle, will depend 
largely on the theoretical accuracy of the ratio £. Note that using |T4jAts| = 0.20 ± 0.02, 
the estimates (J3*7)l and (J3*8)l result into the lower limit ( > 0.81 (at 90% CL.) from 
the BABAR data and the range 0.71 < ( < 1.19 from the BELLE measurement. These 
inferences are not precise enough to distinguish among models of SU(3)-breaking. We hope 
that with the first measurement of R exp [(p, uj)^ / K*j] having been already posted [1], the 
ratio ( will receive a renewed theoretical effort, in particular from the lattice community. 



6. Current and potential impact of R exp [(p, to) 7/^*7] on the CKM unitarity 

triangle. In this part we present the impact of the B — > (p, u) 7 branching ratio on 
the CKM parameters p and fj. For this purpose, it is convenient to rewrite the ratio 
R\h{{p, o>) 7/1^*7] in the form in which the dependence on the CKM-Wolfenstein parame- 
ters p and fj is made explicit: 



i? th [(p,^)7/iT7] 



A 2 C 2 (Mf-m 



2\3 



4 {M 2 B - m\* 



2G(p,fj,e^)+G(p,fj,e^ 



7, m 



(41) 



4 (Ml -ml, j 

Here, the function G(p, fj, e) encodes both the LO and NLO contributions: 

G(p,fj,e) = [1 - (1 - e) pf + (1 - eff + 2 Re [G - pG 1 {e) + (? + fj 2 ) G 2 (e)] , (42) 
and the functions Gi (i — 0, 1, 2) are defined as follows: 

Go = [Ag'-AVrycl **, (43) 
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Table 5: Values (in units of 10~ 2 ) of the real and imaginary parts of the functions G{ 
(i = 0, 1, 2), including the parametric uncertainties, are presented for three values of the 
weak-annihilation parameter Ea- 





Re G 


Im Go 


Re Gi 


Im Gi 


Re G 2 


Im G 2 


+0.30 
+0.03 
-0.03 


4.63 ± 3.89 
4.63 ± 3.89 
4.63 ± 3.89 


-0.48 + 1.50 
-0.48 + 1.50 
-0.48 + 1.50 


3.50 ± 7.90 
10.80 ± 5.96 
12.43 ± 5.69 


6.20 ± 4.09 
9.01 ± 3.28 
9.64 + 3.13 


-1.84 + 5.89 
6.10 + 4.23 
7.86 + 4.06 


3.60 ± 3.38 
9.18 + 3.55 
10.42 ± 3.60 



Table 6: The input parameters used in the CKM-unitarity fits. Their explanation and 
discussion can be found, for example, in Ref. [41]. The parameter r]i is evaluated at the 
scale of MS mass m c {m c ) = 1.30 GeV. 



A 


0.2224 ± 0.002 (fixed) 


\Vcb\ 


(41.2 + 2.1) x 10~ 3 


\v ub \ 


(3.90 + 0.55) x 10" 3 




0.736 ± 0.049 


\ € k\ 


(2.280 + 0.13) x 10~ 3 




(0.503 + 0.006) ps" 1 


m 


1.32 + 0.32 




0.57 + 0.01 


V3 


0.47 + 0.05 


m c (m c ) 


(1.25 + 0.10) GeV 


m t (m t ) 


(168 ± 4) GeV 


B K 


0.86 + 0.15 


fs d \j B Bd 


(215 + 11 + 15^ 3 ) MeV 


Vb 


0.55 + 0.01 




1.14 + 0.03 ±0.02+°;J 3 i°;° 3 


R exp [(p, u) 7 /K*>y] 


0.042 + 0.013 


AM Bs 


> 14.4 ps" 1 at 95% C.L. 








Gi(e) = 2G - [A u 


+ eAW] /Cf cff , 


(44) 




G 2 (e) = G -[(l- 


e) A u + eA^] /cf cfr 


(45) 



Numerical values of the real and imaginary parts of the functions Gi (i = 0, 1, 2), and the 
parametric uncertainties, are given in Table |3| The three rows in this table correspond 
to the decays B — > p ± 7, B®(B®) — > p°7, and B®(B®) — > wy, respectively. It should 
be noted that the function G(p,f],E) (}4~2*|) is related with the dynamical function AR, 
introduced in Ref. [15] to account for the weak-annihilation and NLO corrections, with: 
G(j>,fj,e) = I%(l + AR). 

To undertake the fits of the CKM parameters, we adopt a Bayesian analysis method. 
Systematic and statistical errors are combined in quadrature. We add a contribution to the 
X 2 -function for each of the input parameters presented in Table |H1 Other input quantities 
are taken from their central values given in the PDG review [38]. The lower bound on the 
mass difference AM Bs in the B®—B® system is implemented using the modified % 2 -method 
(as described in the CERN CKM Workshop proceedings [39]), which makes use of the 
amplitude technique [40]. The B s <-> B s oscillation probabilities are modified to have the 
dependence P(B S -> B s ) oc [1 + *4cos(AM Bs £)] and P(B S -> B a ) oc [1 - *4cos(AM Bs t)]. 
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Table 7: The 68% CL. ranges for the CKM-Wolfenstein parameters, Rb = \/ p 2 + r/ 2 , 
R t = - p) 2 + rf , CP-violating phases, AM Bs and R{(p, u) 7/^*7] from the CKM- 
unitarity fits. 



A 


0.2224 


A 


0.79 


-=- 0.86 


P 


0.10 ^ 0.24 


V 


0.32 


-=- 0.40 


Rb 


0.37 0.43 


Rt 


0.83 


-=- 0.98 


sin(2a) 


-0.44 ^ +0.30 


a 


(81 - 


=- 103)° 


sin(2/3) 


0.69 0.78 


P 


(21.9 


-=- 25.5)° 


sin(27) 


0.50 -=- 0.96 


7 


(54 


+ 75)° 


AM Bs 


(16.6 20.3) ps- 1 


R[(p,uj) 1 /K* 1 ] 


(2.3 - 


- 4.3)% 



The contribution to the x 2 -function is then: 



1 1 - A 
Erfc" 1 I - Erfc 



2 y/2 



(46) 



where ^4 and 04 are the world average amplitude and error, respectively. The resulting 
X 2 -function is then minimized over the following parameters: p, fj, A, Bx, r)i, r] 2 , r] 3 , 
m c (m c ), m t (m t ), i] B , f Bd ^B Bd , ^. Further details can be found in Ref. [41]. 

We present the output of the fits in Table d where we show the 68% CL. ranges 
for the CKM parameters A, p and f/, the angles of the unitarity triangle a, (3 and 7, as 
well as sin(20j) with (pi = a, /3, 7, and AM Bg . The allowed profile (at 95% CL.) of the 
unitarity triangle from the resulting fit is shown in Fig. El as shaded region. Here we also 
show the 95% CL. range of the ratio -RgxpKp,^) 7/^*7] = B C x P [B — > (p,u) -y]/B cxp (B — > 
K*ry), which is used as an input in the fits now. We find that the current measurement 
of R exp [(p, to) 7/^*7] is in comfortable agreement with the fits of the CKM unitarity 
triangle resulting from the measurements of the five quantities (Rb, ex, AM Bd , AM Bg , 
and a^K s )- The resulting contour in the p — fj plane practically coincides with the shaded 
region, and hence not shown. We conclude that due to the large experimental error 
on /2 exp [(p, uj) j/K*j], but also due to the significant theoretical errors, the impact of 
the measurement of B — > (p, 07)7 decays on the profile of the CKM unitarity triangle is 
currently small. How this could change in future is illustrated by reducing the current 
experimental error on R exp [(p, u) 7/^*7] by a factor 3, which is a realistic hope for the 
precision on this quantity from the B-factory experiments in a couple of years from now. 
The resulting (95% CL.) contours are shown as dashed-dotted curves, which result in 
reducing the currently allowed p — fj parameter space. This impact will be enhanced if 
the theoretical errors, dominated by A(/(, are also brought under control. 
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Figure 3: Allowed p — fj regions following from the six measurements €k, AMg d , 
AM Bs , a^Ksi an d ^exp[()°) w ) corresponding to 95% C.L., with the dot showing 

the best-fit values. The shaded region shows the current profile. The two outer (solid) 
curves give the 95% CL. constraints in the p — fj plane from the current measurement of 
-Rex P [(p, w) j/K*j]. The inner (dashed-dotted) curves are the 95% CL. constraints from 
an assumed measurement of R exp [(p, uj) y/K*y\ having the current central value but the 
error reduced by a factor 3. The contour shows the potential impact of this assumed 
measurement in the p — fj plane. 

7. Summary We have studied the implication of the first measurement of the aver- 
aged branching fraction B exp [B — > (p, uj) 7] by the BELLE collaboration for the CKM phe- 
nomenology in the SM. Updating the earlier theoretical calculations [15], carried out in the 
QCD factorization framework, in which several input parameters have changed, we have 
calculated the averaged branching ratios for the exclusive B — > (K* , p,u) 7 decays and 
the ratio R t h[(p, 7/^*7]. Using the CKM-Wolfenstein parameters p = 0.17 ± 0.07 and 
fj = 0.36 ±0.04 from the unit arity fits [25], we find B th [B {p,co)j} = (1.38±0.42) x 10~ 6 
and J2th[(p, to) 7/^*7] = (3.3 ± 1.0)%, to be compared with the experimental numbers 



B exp [B -> My] = (1.8lg;|j ± 0.1) X 10" 6 , and R exp [(p, uj) 7/^*7] = (4.2 ± 1.3)%, 



respectively. We see a quantitative agreement between the SM and the BELLE mea- 
surement in the radiative penguin b —>■ d transitions. Leaving the CKM parameters p 
and fj as free, we determine (at 68% CL.) 0.16 < \V td /V ts \ < 0.29 (at 68% CL.), which 
is in agreement with but less precise than the corresponding range from the CKM fits 
l^td/^tsl — 0.20 ± 0.02 [25]. This is, however, expected to change as the experimental 
precision on the branching ratios and R exp [(p,uj) j / K* / ~y] improves. We emphasize that 
the measurement of R exp [(p,uj) y /K* / -f] provides the first direct determination of the ra- 
tio I Vfrf/ Vt s I in rare £>-meson decays. We have also presented updated estimates of a 
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number of the isospin-violating and SU(3)-violating ratios and CP-violating asymmetries 
in the B — > (p, uj) 7 decays. Their measurements will either overconstrain the angle a of 
the unitarity triangle, or they may lead to the discovery of physics beyond the SM in the 
radiative b — *■ d*y decays. 
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